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._ELIMINATION BY INDETERMINATE MULTIPLIERS. 


_ As we have made frequent use of the method of elimination by 

means of indeterminate multipliers, we shall give the panne: of the 
| process, 

The solution of an extensive class of problems may be reduced to 

the determination of —. maximum or minimum of a function U of 


n variables, x,, 2, ...... ; these being in general not independent 
but subject to certain etebileb of condition, 


The method simplest in theory, is to determine, by means of the. 
r equations of condition, r of the quantities r, as for instance, 
x,, in terms of the remaining 2’s; substituting these 
values, U will now be a function of the 2 — r independent vari- 
ables a,-,...... ,; and making the coefficient of each differential 
in dU equal to zero, we get n —Tr equations, which, with the r 
But in practice this method is rarely possible, and still more rarely 
convenient. Instead of this, we may consider each of the n 
quantities # to vary in U, their variations being no longer inde-— 
pendent, but subject to the condition that in their varied state they 
must satisfy the r equations L = 0, or that L, +dL,=0, &c., 
are =(, the variations must be such as to 
: satisfy the r equations 


From these equations we may determine r of the variations in 
terms of the (n — r) remaining; and substituting these in dU, we 
shall have, as before, the coefficients of the ( mn — 7) variations to 
equate to zero. | 
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Elimination by Multipliers. 
The simplest method of effecting this. is the following. The 
quantities d L,, ess as well as dU, are linear in dx,...... 
being of the torm- | | 
| di, aL dL | 
| dx, dx, dz, / | | 
Multiplying them by d,, A,, ...... A, respectively, and adding them | 
to dU, we have . | 
dU + dA, dL, +A, d1,...... = 0, 
an equation which will be of the form | 
+ ...... + Mode. =0, 
in which each quantity M= +r, — + 
we determine the r by that they 
must make the coefficients of the variations 
equal to zero, that is to say, if we determine them from the equa- 
tions M,=0..... M,=0, we shall have effected our purpose. 
For by this means the variations Sap eee dx, are eliminated, and 


there remains 
M,4,d2,,, + M — 


and as these m — r variations are independent, their coefficients” 


must be equal to zero. 3 
Hence, if we determine ,, Ag, ...- A, from the equations 


we shall also have ns | 


Thus the 7 quantities and n x the 
equations 


M, = 0, M, = 0, ‘cain M =0, L, = 0, ...... L, = 0. 


And as it is indifferent which of the variations we suppose to have 


been eliminated, and which set of the - equations M = 0 to have. 


been used for the determination of the \’s, the most general way 
at stating the result is, that we have the » +7 equations 


dU di. ) 
+. +r 2 de, 

dU dL, dL. | 

to determine the n quantities. 2, Ay Ar 
If U be homogeneous in a,...... if L,, L,...... consist 


of terms of not more than two different dimensions, There exists a 


simple relation between the quantities 4, which is very useful in 
many problems. | 
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Suppose U to be homogenedus of dimensions, L,, L,, 


consist of terms of two different dimensions, so that the equations 


L, =0, may be put under the form 


M,+A=0 N,+B=—=0, &c. 


where M, is homogeneous of a dimensions, and A a constant. 


Then multiplying the equations (1) by 2, a, ...... and adding, | 
we have, by the property of homogeneous functions, es 


) 
0 = uU + aM, + + &e. 


or uU = aA), + OBA, + cCh, + &e. 


We shall now proceed to ib a few examples of the application 
of this method. 

If we take the tangent plane at any point in a surface for the : 
plane of ay, and the normal for the axis of z; and if.7, s, t, accord-— 


ot. 

dx*’ dx dy’ dy?’ 
at the point in question, a, } the cosines of the angles that a plane 
passing through the axis of z makes with the planes of «xz, yz, and 


ing to the usual notation, represent the values of 


if p be the radius of curvature of the section of the surface made 


by this plane, it is known that — 


with the relation = Lo (2). 


If we wish to find the direction of the section of greatest or least 
curvature, we we substitute YW] — a’ for Bs we have then 


+6 


from which we eee determine the maximum and minimum values : 


of p with the corresponding directions. 


But without eliminating 3, we may consider it as a function of a 
given by equation (2), and make a and # vary in (1), subjecting 
them to the condition of satisfying ( an in their varied state. 


Hence 
+ $3) da + (sa + (3), 
If we eliminate dj}, the coefficient of da equated to zero will 
give a relation between a and (2; to do this, multiply (4) by A, and 


add to (3) 


(ra + 83 + a) da + (sa + {3 + AB) dB = 0. 
If, then, we determine X from the condition 
Sa +AG=— 


we must also have 
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52 Elimination by Indeterminate Multipliers. 
These equations give us at once 


| which shows us, that the two directions of maximum and minimum 
curvature are at right angles to each other. 
Multiplying by a and B, and adding, 
4 + = + B’) 0, 


Substituting this ae (5) and (6) may be put in the form 


(- - ) 
(~-r) 46, 


and multiplying together, 


The relation --+ A =0, corresponds to the relation uU =aar, 
| + &c., which we have found in the general case. 


It may be observed in this case, that if we merely wish to deter- 
mine the relation between a and 6, we may eliminate dj at once 
and obtain the relation; and if we have three equations and three 
arbitrary variations, we should get tiie relation by eliminating two 


of them by the method of cross multiplication. An example of this 


will be found in the investigation in p. 7 of our last Number. . 


For a second example, let us investigate the properties of conju- 
gate diameters in surfaces of the second order. | 


Let | 
x* y ge | 


be the equation of the surface referred to any three conjugate di- 
ameters, whose lengths are a’, b, c. 


If J, m, m be the cosines of the angles which the 8 axes make with 
each wa r any radius, 


= + 4 4 Qyz + + Qnyz......... (2). 


When + coincides with one of the principal axes, it is a maximum 
or minimum, and we have the condition dr = 0. 
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If we determine z from (1), and substitute in (2), r is a func- 


tion of x and y only, and hence we should have 


| which, with (1) and (2), will serve to 


This would be very difficult ; but we may, as DerOres differentiate | 


(1) and (2), L,Y, z to- 


If we multiply he second by ’ and add to the first, and js 
mine A from the condition 


tly =o, 


we have an aoe in which dz has been eliminated ; and as 


~ dx, dy are independent, their coefficients must also be equal to 
zero, OF 


nz =0, 


If we multiply these by 2, y, z respectively and add, we obtain, - 
in consequence of (1) and (2), 
= 0, 
Substituting this value of A in the equations, they become 


mx 4+ ly + (1-5)z=0 


Eliminating x, y, z between these by the method of cross s multi- 
plication, the resulting equation is 


+ = 0, 


a cubic equation in 7*, the roots of which are the squares of the 
principal axes; and by a discussion of which the properties of the 


conjugate diameters may be deduced.—Vide Hymers’ Analytical 
Geometry. 
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54: Linear Equations of Finite and Mixed Differences. 
As another example, find the values of x, y, 2 which may make. 
the function 
(x 1) 1) a maximum, 
Ys being connected the equation, | 


Taking the gers differentials of both equations; 
dx dy dz | 


2 
0 = dx log a + dy log b 4 dz log c. 


Multiply the second by —A, add and equate to 0 the coefficients of 


- each differential, then 


] 

A log 6b = 
Multiply by xe+tl,y Zz and add ; then 
A(wlog a + ylog + loge + log a + log + log e) = 3, 
or X log (Aabe) = 3. 


5 log c= 


log a = 


Therefore substituting 
log (Aabc) 


3 log a 
therefore. 


Aabe) 


$loge 


log (Aabc)? 
log (abe)> 
M.S. 


Il._ON THE SOLUTION OF LINEAR EQUATIONS OF 
FINITE AND MIXED DIFFERENCES. 


In the preceding Number the principle of the separation of the : 


symbols of operation from those of quantity, was applied to the 
solution of Differential Equations. We shall here apply the same 
principle to equations of Finite and Mixed Differences; but, as the 
method is not very different from that previously delivered, there 
will be no necessity for dwelling long on it. ~ 

The general form of the linear equation of Finite Differences 
with constant coefficients, is 


+ + + Ru,,, + Su, X, 


where A, B, C, &c. are constants, and X a function of a By a 
known relation we can express the quantities u,, v,,,_,, &c. in 


q 
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terms of u, and its successive differences, so that the equation may 


be transformed into one which may be i under the form 


which we might proceed to solve by means of the separation of the | 
sy mbols. But it will be more convenient to proceed in the follow- 


manner. Since we generally, 


= (1 4 


we » may consider 1+ Pye as a separate symbol, subject to the same 


~ laws as d and A, the effect of which, when applied on 1,, is to con- 


vert w, into v,,,. If we represent, for the sake of shortness, 1 Ba he A 


by D, ‘the wi sone will take the form 


+ AD™! 4+ &. + RD + _X, 
or, as we may write it for convenience, 

| F(D) u, = X. 
Then, converting F(D) into factors of the form D — a, a being 
one of the roots of the equation F(z) = 0, we. can perform the 


inverse operation for each factor separately, and so arrive at the 


value of wx. To do this more readily, we shall avail ourselves of a 
theorem similar to that given in page 25. This theorem i is, that 


(D - = antran (Xa-*), 


For if we expand the first side, it may be put under the form 


| 
— &e.) X. 


Now DPa-z — and a~P =: 


so that the expression may be put under the form 


(D"D‘" 4 &e. ) Xa~*, 


where the accented letters refer to q-2 and the unaccented to X. 


This expression is equivalent to 


— (A A’ + AA’ yr -= (Xa™*), 
by a theorem analogous to that of Leibnitz. 


To apply this to the general equation F (D)u; = X, suppose 
F (D} to be resolved into binemial factors, such as D — a, so that 
it becomes 


(D — (D By) (D — Gy) Uy = 
Then taking the inverse operation of D — a, ans 
(D) u, = (D 3 
representing the product of the n — 1 binomial factors by I’ (D), 


and omitting the arbitrary constant, as included in the sign of in- 
tegration. Again, performing the inverse operation of D —a,, 


we have 


k 
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\\a, 
and this by parts, 


xr—1 


a 
Proceeding i in the same manner with the next factor, 
(a, — (a,—a,) (a, — a) (4, — 45) 
+ * 
(a, — dy) ( 


after integration by parts and reduction. And in the same way 
we may proceed till all the binomial factors of operation are ex- 


hausted. It is obvious from this, how close an analogy this result 


bears with the corresponding one in Differential Equations. 

If there be x equal roots, the theorem given above enables us to 
arrive meee at the solution. For in this case we should have 

or constants, which were before neg- 
lected, | 


Uy = at * =" (Xa "(C+C2x +. Ca"). 
_ Linear equations of the first degree, with variable coefficients, may 
be easily converted into equations with constant coefficients, by a 


change analogous to the change of the independent variable in the 
Differential Caleulus. For let 


Urs, + aP = 


be the given equation, where P,; is a function of x We may 
assume P, = ——, so that dividing the becomes 


z 


which may be solved as a linear equation in- , from which, om 


Q,’ 


found, we may determine uw, The form of Q,,, is evidently 


The same method may be sometimes applied to equations of 2. 


higher order. If, for instance, we have the equation 
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| Dividing. by ee we get 


which may be solved as an equation in 


u, determined. 
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us represent the continued product ¢ (1). 9) + 1) 
by P,,, Then 


(@) ple 41) = and 41) = Pas, 


P, 


—*_, and so the value of 


It will be seen that the method here laid down of itneatinar 
equations of differences, possesses an advantage over the common 


one of not requiring any assumption of a form of solution; but — 


this advantage is even more displayed in solving equations of 
mixed differences. In these the two symbols d and A are in- 
volved at the same time; but as they are the characteristics of 
independent operations, the same principles which were applied to 


each separately will hold when they are combined. Without dis- — 
cussing the general equations of mixed differences, we shall pro- 


ceed to a few of the more interesting examples. Let us have the 


equation 


d dy 
bay + aly | 
By separating the symbols this may be put under the form - 
d | 
(— +) (A4ay=X 


We shall obtain a different result according as we operate first 


with the one or the other. factor ; _ taking then the differential factor 


first, we have 
(Atay=e + Ce 
Now At+a=1+A-(1-a)= 


Therefore, performing the inverse — by the theorem given 
above, page 55, we obtain 


+ $4(1—a)*¢(sin 2ra2, cos 2rx), 


the last term being the ot function arising from the 


operation $D — (1 — 


And performing the indicated in the second term, 
we have 


+ + (1 — ¢(sin cos 


If we begin with the other factor we have 
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d 
and y =e fdr (1 — a)? 
| + [dx (1 — a)? + Ce, 
_ where ¢ is put for @ (sin 27a, cos 272). 


The next example we shall take is one given by Paoli, where 
two variables are involved, 


d 


which may be put under the form 


| where the D refers to the z only. 
Then, following the same principle as that pursued by Mr. Great- 
heed in the Philosophical Magazine for September 1837, that is, 


considering the characteristic rs as an independent quantity, we 


an arbitrary function of y being substituted for a constant, as the 
operations are partial. This may evidently be put under the form | 


dy* dy* 


If P,. y =O, the result is the same as that given by Sir John 
Herschel, in p. 38 of ‘his Examples. He likewise gives the 


equation | 
(2) b =0 
which may be put under the form | | 
d? 
D? — aD — b ) 
| ( dy af; 

where D affects x only. Or, putting it into the form of factors, 


where m, n are the roots of the equation 


d d 
Then (D n =) u,y=m =) GY), 


have 


Us, y 
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as the functions of y are quite 


In the same manner we might integrate this equation, if the 
second side were some function of 2; and also equations of the 
form 


d 

— 
where A affects z only. But it is needless to dwell on ‘cien and — 
we shall therefore proceed to the integration of an equation, in 
which the advantages of this method are very conspicuous. Mr. 
Airy, on the hypothesis that the distances of the particles of the 
luminiferous ether are not infinitely small compared with their 
sphere of action, has given the following ones for — 
the disturbance of a particle : 

at 

de * = A ho? 
where A affects x only. This may be put under the form 


Following the same method as before, by splitting the operating 
factor into two, and integrating with each separately, we obtain 
A 
where (x) and are arbitrary functions of 2. To reduce 


this to a more intelligible form, we may avail ourselves of Fourier’s 
formula for the transformation of functions. From it we have - 


(2°) = f_ dag(a) dp cos p(x — a). 
Substituting this form, 


A?u 


a 
| da g(a) dp cos p (a — a) 
A 
dp ‘Vita cos p (a — a). 
We must now expand the exponential, and operate with each 
term separately on the cosine. When expanded it becomes 


a 
(1 + + + Xe.) cos p (4 — a). 
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60. ~=Linear Equations of Finite and Mixed Differences. 
Now | 


therefore Vita cos p (@ — - a) — 2sin sin p(x — a), 


sp (x )= 2 ( 
a -( sin “> 


and so on for _— repetition of the operation, so that the series — 
becomes | | 


2 22 4 | ‘ 
_ (sin) (2 sin?) — 


9) * 781.234 2 


X cos p(«x# — a) 


(28 sin 9 193 2 sin 4. se, p (2 a) 


t 
= cos (=" sin? ) COS 4) = sin sin sin p a). 


Similarly for the other term we shall obtain 2 : 
Qat . ph  (2at Or 
cos sin cos p(x — a) + sin (> sin sin p(« — a); 
and if these expressions be substituted in the expression for w;, 4 
we shall obtain a result free from the symbol A: but a slight 


- change in the form of the function will reduce it to a somewhat 
better form; for if we make | 


the expression becomes 


da f(a) dp cos sin cos p (x — a) 


da F (a) dp sin (= sin =) sin p (a — a), 


which is the general solution . “ ih equation. The form of the | 
arbitrary functions are easily determined from the initial cireum- 
stances. For if in the last expression, or in the original integral, 

we make ¢ =0, we find 


(x), 
which determines the form of the function f from the initial state 


of the particles. If we differentiate the expression for u,,; with 
regard to ¢, and then make ¢ = 0, there results 


dt = = E (a4 


which determines the form of the function F from the initial 
velocity. 
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Mr. Airy has found a particular integral of this — by as- 
| suming as a form of solution 


and then. determining v. But his result may be arrived at, merely — 


by assuming the form of the arbitrary function to be that used by 
_ Taylor in the solution of the problem of vibrating chords, and 


which is usually taken as the form of the function in the undu- — 
latory theory, and is for that reason well —_— for comparing — 


the results. Let then 
= — at), == 0, 


[effecting the operations, in the original intemral by the 
same method as that used in reducing Fourier's function, we 


obtain | 
A 


This expression gives us for the velocity of propagation of the 
wave, which is the coefficient of ¢, 


which is no longer independent of ), and therefore the velocity will 


be different for different values of X; and as the index of refraction | 


is inversely as the velocity, it will be nearly inversely as A, and 
will thus be greater for violet than red light. This result then, as 
far as it goes, accounts for the phenomena of dispersion. 


IIT. —DEMONSTRATIONS OF SOME PROPERTIES OF 
THE CONIC SECTIONS. 


1. Tue position of the circle of curvature at any point of a 
Conic Section, may be readily determined by the following con- 
struction. Describe a circle touching the curve at the given point 
and cutting it in two others, then the chord in the conic section 
which passes through the given point and is parallel to the line 


joining the two points of section, is a chord of the circle of curva-- 


ture at the given point. 


Taking the origin at the point in the curve, the equation to the 
conic section is 


y’ + Bay 4+ Cx? + Dy + Ex = 0. 


4 
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If we make the normal the axis of x and the tangent that of y; 
when x =0 the two values of y must also =0, and therefore 
D = 0, which reduces the equation to | 
| + Bry 4 Co? + Ex = 
re equation to a circle touching the curve at the given point, : 
‘that is, being also a go to the axis of y at that point, is 
2? + —.0. 
At the intersection sf the curves we may combine the equations 
in any manner; subtracting them we have. : 


E') = 0," 
which gives 
‘The first of these equations gives the origin, the second gives a. 
relation between the coordinates of each of the points of section ; 


and as it is linear it is the equation to the straight line joining 3 
them. Now the angle which this line makes with the axis of a is 


determined by the ratio 7 : ; and as this is independent of the 


circle, the line will remain parallel to one position, as the circle 

varies in size. But when the circle becomes the circle of curva- 
ture, one of the points of intersection coincides with the point of 

_ contact, and the line joining the points of intersection, which re- 
mains parallel to one position, must pass through the point of con- 
tact. Therefore, if from this point a line be drawn parallel to a 
known position of the line of intersection, it will pass through the 
point in which the circle of curvature cuts the curve, and will thus 
be a chord both in the curve itself and in the circle of curvature. 
Knowing now one chord of the circle of curvature, and the posi-. 
tion of its diameter, which coincides with that of the normal, we 
can determine the circle altogether. 


2. The following is another and a very curious ssaibini of deter- 
mining the centre of curvature in the Conic Sections. It was first 
given by Keill, but seems to have been rather strangely neglected 
by the subsequent writers on this subject. 

From any point P (fig. 1) in the curve draw the normal PN 
cutting the axis in N; at the point N draw NQ perpendicular to 
the normal, and meeting the focal chord through P in Q. From Q 
draw QO perpendicular to the focal chord and meeting the normal 
in O; then QO is the centre of the circle of- curvature at the point 
P. Draw FY from the focus perpendicular to the tangent, and let 


In the right-angled triangle PQO, we have plainly 
PO.PN. = PQ’. 
But from the similar triangles QPN, FPY, 
PQ.FY = FP.PN ; 


ey 
- 
2 
| & 
€ 
‘ 


Demonstrations of some Properties of the Conic Sections. 63 


therefore PQ? = N? , and PO=N 


_-Also in a conic section we have, if R be the radius of curvature, 


3 3 

mr 

R = = also —~, 

where mis half the latus rectum. 
From which we have — = -—,» and therefore 
m | 


m* 


and is contre of the circle of curvature at 


3. In the Cambridge Transactions, Vol. 3, Mr. Morton has de- | 


monstrated a number of curious properties of the Conic Sections 
in relation to the generating cone; but -he does not seem to have 


noticed the following one. If a sphere be described round the 


vertex of a cone as centre, the /atera recta of all sections of the 
cone made by planes touching the sphere are equal. Taking the 


vertex of the cone as the origin, and the axis of the cone as the 
axis of a, the equation to the cone will be 


And if we suppose the cutting ae to be einer to the 
plane of az, its equation will be 


where r is the radius of the sphere, and a the angle which the per- 


pendicular from the origin on the plane makes with the axis of z. 
Eliminating z between the two equations, we get 
a? (cos? a — m sin? a) + y? cos* a + Qm?rex sin a = m?r?, 
which is the equation to the projection of the section on the plane 
of xy. This equation, which is that to an ellipse, is not referred to. 
its centre, but if we so refer it, it becomes | 
a? (cos? a—-m? sin? a)? +. y” cos? a (cos? a— m? sin? a) = m*r 2 cos” a a. 
Now, if a’, b' be the axes of the projection, 
mr? a m*r? 
(cos? a — m? sin? a)? ~~ €0s? a -— m? sin? a 


If a, b be the axes of the section, as the cutting plane is perpen- 
dicular to the plane of wz, and makes an angle a with the plane of 


XY; 


a@d=acosa 
And for the ldatus rectum, 
Qb2 


a 
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~ which being: independent of a, is the same for all sections for which 
ris the same; that is, for all those which are made by planes touch- 
the sphere. | | 
From this it appears, that the datus rectum is equal to the di. | 
ameter of the sphere multiplied id the —_—: of half the vertical 
~ angle of the cone. 


-D. 


THE MOON'S. ANNUAL AND SECULAR 
INEQUALITIES. 


THE following is a simple method of obtaining the expressions for 
the Moon’s annual and secular equations, and which may be applied 
to determine those parts of the Moon’s inequalities in longitude, 
which depend on the central disturbing force, and are independent 
of the eccentricity of the Moon’s orbit and the square of the dis- 
turbing force. 

In this investigation the Moon is supposed to be moving ina 
circle, the radius of which is variable. The equable description of 
areas is independent of the magnitude of the central force. Hence 


h = (pa)? is constant, p» and a varying. 


Let cu, ca represent the finite variations of p and a at any time. 
We suppose their squares to be neglected in this investigation, and 
hence they may be treated as if they were infinitesimal. 


Differentiating the logarithm of the expression for hy 


aT 6a 
5s 
6d6 op 83. 


and substituting the a. of from equation (1), and putting 
dcé for (éd6), 


fe 


which will give the Moon’s inequality in longitude. 
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~The central force is by hypothesis + ; in n the case of 


Moon disturbed by the Sun, this is equal to f. ~—i “4 , where 


_m’', r’ are the Sun’s mass and distance. 
‘a? 

As the Sun describas an elliptic ath about the Earth, 
+ € cos —w') 
a(l-—eé*) ° 
Let m be the ratio of mean motions-of the Sun and Moon, (3 the 
_ Sun’s longitude when the Moon’s longitude is 0, 


0' = + + terms multiplied by e. 
ig Hence, as far as the first power of é, 


a 


r 


substituting this value, and putting m* for —.—, | 


cos (m0 + (2 ). 


The effect of the first term of this.expréssion is to increase the 
Moon's periodic time by a constant quantity, which is in fact 
included in the mean motion determined from observation. The 
second term will give the inequality depending on the place of the 
Sun and the eccentricity of its orbit. Therefore, confining 0 to 
this part, 


— {cos (md + B — db, 
| — 3me sin (md + — w), 
= — $me' sin (Sun’s mean anomaly ); 


which is the expression for the annual equation. Vide Airy, 
Lunar Theory, No. 59. 


The effect of that part of the Sun’s attraction which is inde- 
pendent of the Moon’s place in her orbit, is to diminish the Earth’ 5 
a... attraction by the quantity 
| ma 


This depends on the Sun’ S position in its orbit. 


Its mean value will be found by adding its values taken at every 
instant during a complete revolution of the Sun, i.e. a year, and 
dividing by the number of instants ; this amounts to integrating the 
quantity, multiplied by the element of BME, for a year, and di- 
viding by the length of a year. 
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Let T represent the length-of a year. 


ma 

Ona”? m aft 
Hence mean value 
2 Q 
(1 - 


Neither the Sun’s mean distance, nor the length of the year 
which depends on it, alter sensibly. But in consequence of the 


action of the planets, the eccentricity of its orbit e¢ has been con-— 


stantly diminishing for centuries. The expression we have Just 


found shows, that in consequence of this the value of the Sun's dis- _ 


turbing force is constantly diminishing, or the Earth's attraction 
increasing ; and hence arises that inequality celebrated in the his- 


tory of Astronomy as the secular acceleration of the Moon’s mean 


motion. 


To find its amount, suppose E the eccentricity of the Sun’s 
orbit at any epoch, 


The effect of the first two terms is included in the Moon’s mean 
motion, as determined by observation at the epoch in question ; the. 
inequality in longitude at any other time, depending on the dimi- 
nution of é, will be 


o= =m? e'2) dé, 


T 
| 
3 ‘9 
— m*— _me* + &e. 
| 
; 
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which is the equation given in o Mécanique Céleste, liv. vii. 
No. 15. 


| Let E’ be the annual diminution Py é’, and let the year be taken 
as the unit of time; 


eo E — Et 4 — QEE't nearly, 
md@ = 2rdt; therefore 
60 = 3m7EE’ J 2tdt, 
= 3mrEE'é, 
3mr(ecc’. of Sun's orb. )(ann.v var. of ecc’. )(sq. of n°. of yTs. ). 


H. T. 


V.—ON THE EXPANSION OF A FUNCTION OF A 
BINOMIAL. 


THouGu Taylor’s theorem natively presents itself at the very 
entrance to the Differential Calculus, many points relating to it 
‘are imperfectly understood even by the best. mathematicians. It 
does not seem to have occurred to any one, that the particular 
expansions of f(a +) in fractional or negative powers of h, to 
which we are obliged to have recourse where Taylor's series fails — 
when a2 becomes a, are derivable from another general expansion: 
of f(2 +h). This we shall show to be the case. It seems, 
indeed, that many persons, misled by Lagrange’s proposition, “ the 
expansion of f(a +h) cannot contain any fractional or negative 
powers of A unless a have a particular value,” think that there is 
no such other general expansion; but the above proposition ought — 
to be enounced thus, “that expansion of f(a +h) whose first 
term is f(x), cannot contain any fractional or negative powers | 
of h, and cannot fail unless x have a particular value ;” for the 
propositions, “the first term of the expansion is f(2),” and “the : 
expansion cannot contain fractional or negative powers of fh,” are — 
reciprocally dependent on one another, as will be seen. Professor 
Peacock, in his Report referred to in page 11 of: our last Number, 
_ has employed general differentiation to express the general term 
of the expansion of f(a +); but in all other respects we believe 
that the contents of this Article are new. 


2. Let the following general series be assumed for f(a +A), 


where the indices a, a’, a’ are supposed to be in order of magni- 
tude, either increasing or decreasing. Equating the first differen- 
tial coefficients with respect to x and hf, 
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| de dx dx 
aAhe-) 4 4 ANN | 


One way by which this equation may be satisfied is, by making 
a =); this will lead to the known series of Taylor. But there is 
another way by which equation (1) may be satisfied; this is, by 


dA” 


~~ making _ equal to zero. Others of the differential coefficients 


may or may not vanish; we will suppose first, for the sake of: sim- 
_ plicity, that they do not. Comparing the other terms with respect _ 
both to the indices and coefficients, : ; 


a =a — a’ =a — 1, &e. 


and since A is a constant = Cy. 
| 


a (a — 1) a—l 


Hence 


f(a +h) = Coe (< Cyr + 


AS ( ( ) 4 4 4. eee (2). 


3. It is only in a particular form of f (x) that f(a +h) can be 
expanded in this series. To determine that form, put x = 0 in the 
last equation, whence 


and changing A into a, 


f(a) = Cyt 4 4 

It follows then, that when f(z) is of, or can be expanded in the 
form (3), which includes all the simplest functions for which Tay- 
lor's series can fail, there exists another general development of 
f(x +h) besides Taylor’s series. In fact, the series (2) is the 
same as would result from changing x into x +A in each term of 
(3), and developing each term by the binomial theorem in descend- 
ing powers of A. It is evident that the series (2) cannot fail for 
any value of x. If all the indices of x in (3) be positive integers, 
the series (2) will be identical with Taylor’s (which will in this 

case terminate) written in the contrary order. 


4. The coefficients of the several powers of h in (2) may be ex- 
pressed in terms of the general differential coefficients of f(x) with 
respect to a. Putting f(#) =, and using the formula (P) in 
page 20 of our last Number, we evidently have from (3), 


a5 
> 04 
| 0 
A" = /(a—1) 
1.2 0 ] 2’ 
5 
| 
4 
| 
2 
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‘As was proved in the article on General Differentiation, . the 


terms in the preceding differential coefficient which contain nega- 
tive powers of a (if there be any such) are finite, but, unless 


a — r be a positive integer, P(a — r) is infinite; hence in dividing 


by P(a — 7) the terms involving negative gs peg of x disappear, 
and we have 


2) —7 + 1) 

which is the coefticient of in equation (2); 
hence when f(x) is of the form (3), f(a +A) may be — 


in the series 
P(a) da P(a —1) P(a — 2) dx*~? ( 
When a is a positive integer, the terms involving negative 
powers of x would not disappear naturally from those differential 


coefficients whose index is positive, but they must be rejected, in 
order that the expansion may agree with (2). 


5. The Series examples may serve to explain the subject. 


bet: ss 7 Va a — w: the expansion of this in descending powers 


of is 


which is of. the form (3), and we find in this case a = 5° Also 


from this expansion we obtain 


| | 2.4.P ( — 2) J 
| 3\ 
leu ] l.] ) 
— f on | i P (5) a- * (? 
dv? 2.4P(— 1) | 


70 On the Expansion of a Function of a Binomial. 
‘Substituting these values in the series (4), obscrving the rel ation 


P(n) = tn — 1), and. that P( —1) = we find 


the (p + 1)!" term 


axP7) 


when x receives the particular value a, the terms after the second 
vanish, and the series becomes 


] 
—1( + ah’). 


6. Again, let u or f(x) = /a?— a Since this is equal to 


22 
we have a= I], 
du. 
je 
dz? 234 201.2 
and since the values of P(1), P(0), P(—1), P(—2), P(— 3), 
at 91.0 respectively, we have, rejecting negative 
powers of x from and w, 


the (p +1)!” term being (provided p be not less than 2) 


bel (p— 2) (p— 8) 

aizp-4 
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This will agree with the expansion of 2 + Yah, when x has the 
value +a. | 
7. Returning to equation (1), let us suppose that several of the 


terms in its first member, besides ' vanish, and that the indices — 


of A in these terms are (, y, &c. By comparing the other terms 
with respect to their indices and coefficients, we should obtain for 
the expression of f(a + h) the sum of the series in (2), and of 
others similar to it with (3, y, &c. instead of a and different con- 
stants. Also, if we denote the second member of (3), without re- 
gard to the values of the constants, by @(a, a), we should find 
that the form of f(2), corresponding to the more general form of 

p (@ +h), is 

g(x, a) + B) + &e.; 

d which form includes all functions that are commonly called Algebraic. 

We shall always suppose the difference of any two of the indices 

Bp, y, &e. to be fractions; for if it was an integer, the two series 

would unite into one. 


8. The more general equation to (4) will be 
P(3) dx" P(g—1) ae 


&e. 
with this provision, ‘that whenever the index of differentiation is a 


positive integer, we reject from the differential coefticient terms 
involving fractional powers of 2. 


For since $B — a is a fraction, the index of x in every term of 


a » (x, 3) is fractional, and therefore its coefticient finite, so that 
- | a = terms will disappear, together with the negative powers of 
r in 5 Ta? (x, a), on being divided by P (a), which is infinite when 


a is oS a positive integer. Hence 
| 
| dat da ? (2, a), 
and for the same reasons en 
dx" ? (x, a), dx (x, 3), &e.; 
so that, after expanding each of the functions 


o(x+h, a), 3), &e. 


we obtain equation (5). 
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9. As a function which requires the more general expansion (5) 
we may take (x? — a*)*. This expanded in descending powers of 
x by the binomial theorem is 
$55, Sl 5-3 
x ax + au x +.. | 
, whence we see that the quantities a, i Bs y, &e. are three i in sinha, 
and that their values are 1, 4, —4, for the rest of the indices in the 
last series differ from these by integers. The rest of the process is 
similar to that used in § 5 and § 6. 


10. We have not yet obtained the most general expansion for 


+h). Suppose that f(x) = Vx —a, and that we wish to 


expand f(a +) ina series that shall not fail when «=a. We 
may expand the factors e**" and Wx + 4 — a separately, and mul- 
tiply the series together. Now e*t” can only be expanded in as- 
cending powers of h, and the other factor must be expanded in 
descending powers of 4. By the multiplication of the two series, 


: therefore, we should obtain one unlimited both ways, as long as # 
_is not equal to a. | 


11. ‘In order to obtain a . general formula to suit such cases, as- 
sume for f(x +h) a series unlimited both ways, namely, 
+ Bh? 4+ Ah* 4+ Aht + BAS 4 Ch’ + ...... 


where we suppose the indices to increase towards the right, and a 
to be the least positive index. 


Equating the first differential coefficients sige to x and A, 


dB dA 4 aA 
a 
eee + ] h + dx h dx 


+ aAh! + BBAP + ...... 


If we suppose at first, for the sake of simplicity, that only one of. 


the indices in the assumed series lies between 0 and 1,  — 1 must 
be the least positive index in the second member of the last — 
tion, and by comparing the terms we obtain 


6—-a—l, a=a—l, a=fp—1, B=y—1, 


dB dA dA dB 
and therefore, ee the other quantities in terms of A and a, 
f(x + h) = = 1) [2Adx? + [Adx 
a dx (a 1) (a + 2) (6). 


12. It is easy to see what would be the consequence of sup- 
posing several of the indices in the assumed series to lie between 
O and 1, namely, that we should obtain for f(a + A) the sum of 
several series similar to that in (6), but in which A would be 
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replaced by B, C, &c. successively, and a by f, y, &c. Such, 
then, is the most general expansion for f(x +h). It shews that — 


f(a +h), in its most general form, is the sum of several simpler 
_ functions of x -+ 4, which may separately be expanded by the for- 


mula (6). It is enough to consider one of these functions at a 


time ; we therefore return to that equation. 


. Let us examine the circumstances under which the series in 
6)a can become limited in either direction. The form of the series 
will not be altered if we suppose 4*A to represent no longer the 


term with the least positive index, but the first or last term that | 


does not vanish. Now all the terms to the left of that may be 

made to vanish by supposing a = 0, but by no other supposition 
while the value of w is general. This leads us immediately to 
Taylor's series, which we thus see to be the only general expansion 
of f(x +h) in ascending powers of / exclusively ; and from there 


being no quantity in it left undetermined, such as a in (4), to be 


applicable to all functions. 


The terms on the right of A*A may be made to vanish by sup- 
-posing A a constant, and by that supposition only. The equation — 


then becomes the same as (2). 


14. It remains to determine the value of the quantity A in — 


equation (6). This may be done in the following manner : 


ee y be a function of x of the form (3), z any other function, 
-and let y' and 2 represent the values of these : quandihes when x 
x-+h. Then by (4) 


Y= P(a—1) * P(a—%) daw? 
and, by Taylor’s series, | 


h2 
l dr 1.2 dx 


The multiplication together of these two series will produce one 
unlimited in both directions, in which the coefficient of h* will be 


P(a) P(a—1) dx * 1.2.P(a—2) aaa 
dry a a(a—1) d**y d*z | 
* 1 da 12° * aa * 
which, by Leibnitz’s theorem generalized, is equal to — 
P(a) 
In the same manner it may be seen that the coefficient of h*~? is 
d*—Pyz 
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and that of | | 

P(a + p) daxtP” 
if it be remembered that all the differential coefficients of y whose 
indices are greater than a are finite. The terms of the series, 
therefore, follow the same law as those of the series in (6); and if 


we suppose yz = u, the two series must be identical, so that by de- 
termining the quantity A, (6) becomes 


f(x +h) 
P(a—1) P(a) da P(a fl) ( 


When it requires the sum of several series such as that in (6) to 


express the value of f(a +h), the determination of A and the 
other corresponding quantities presents some difficulties, so that we 
must ieave it to a future occasion, and perhaps to-an abler analyst. 


_VL—ON SOME ELEMENTARY PRINCIPLES IN THE 
APPLICATION OF ALGEBRA TO GEOMETRY.* 


In pure Algebra, the independent existence of signs of affection is 


the immediate consequence of the universal applicability of its — 


rules of operation. Thus, if the operation indicated in a — 6 is 
always possible, whatever the relation between a and 6 may be, it 
follows that a (a+c), a—a-—e, or (since a—ad= 0) 
must be a possible result of such an operation. The negative sign 
of affection, therefore, in Algebra, although it represents something 
not implied in the primary arithmetical definition of subtraction, is 


yet intimately connected with that operation, and can have no > 


meaning distinct from that of subtraction in any case where sub- 
traction becomes possible, or, in other words, a +(— 6) and 
a — (+ 5) are identically equivalent. 

But in applied Algebra the case is somewhat different. There 
may be quantities of specific natures, which by their reciprocal 
relations afford us an obvious interpretation of the signs of affection, 
and so enable us to define Addition and Subtraction as operations 
‘inamanner no less general than is done in Algebra itself; that is, 
we may lay down rules for the performing those operations, which 
rules shall not be limited, as in arithmetic, by the particular rela- 
tions of the quantities, but shall be applicable in all cases. In such 
cases it is clear, that the independent sign of affection is the found- 


* From a Correspondent. 
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ation, and not, as in Algebra, the consequence of the rule of opera- 
tion; it exists in the very nature of the quantity which is to be 
subject to the rule, and is the condition which makes such a rule 
possiblew. In such cases, therefore, we shall have two distinct mean- 
ings of the negative sign —, as it indicates a quantity of a particular 
_ kind, or an operation to be performed ; the only condition connect- 
ing the two being, that the adoption of either interpretation must 
lead to such results as may be considered algebraical equivalents. 
There may exist any specific diversity, as in the case of geometrical » 
parallelism ; but such diversity, being of a kind not comprehended 
in the algebraic symbols, forms no obstacle to their algebraical 
equivalence. 
_ In the application of Algebra to Geometry, it is to be borne in 
mind, that the subject represented by the algebraical symbols is not 
- geometrical extension, but rather extension combined with direc- 
tion. It is not the distance between two points, but the distance of 
one point rightwards or —— of another. It is, in short, theo- 
retical motion. — 

Quantities of the former, pure geometrical kind, have no opposite 
relation of plus or minus to their two extremities: subtraction ap- | 
plied to such quantities, is the expunging of a line, and when the 
magnitude has once been entirely expunged, the operation can be 
- carried on no longer. The negative sign cannot be applied to dis- 
tance alone, but to distance or progress in a given direction. A 
man who travels north, and afterwards travels the same or a greater 
distance south, has not traversed no space, nor a space less than 
zero, but his progress north is upon the whole either zero or nega- 
tive. We say, therefore, that the quantities we have selected for 
the application of algebraical reasoning are essentially positive and 

negative, independently of any rules of addition and subtraction ; 
they are quantities which are not adequately defined until we know 
whether they are preceded by + or —, and we now proceed to in- 
quire for necessary rules for the addition and subtraction of quan- 
tities of this description. 

~On Addition. 

Let AB be a positive line drawn from A towards B, and repre- 
sented by a or +a: required to add to it a positive line equal to 
+ b. 

Now all positive lines must be drawn in the same direction, say 
towards the right of the point of starting, or that point from which 
their description begins. To increase AB by a line drawn from 
left to right is impossible, except we start from the point B, or the 
_ final point in the motion by which AB was described. We must 
therefore continue the motion AB in the same direction, making | 
BC = 6, and AC, measured from A to C, and not from C to A, is 
the sum of the positive progressions AB + BC, or a + 6. 

_ If we wish to add a negative quantity, or — 6, to AB, it is clear 
that the operation cannot be performed at A, the origin of AB; for 
a neggtive line drawn through A will increase the whole AB, and 
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if such a line were drawn equal to a in length, the result a — a, 


instead of being zero, would be a geometrical line equal twice AB © 
in length; to add — 6 to d, therefore, we must proceed in the | 


same way as when we add + 5 to a, i.e. we must draw BC from B 
in the negative direction, and as before the sum of the two is AC, 


a geometrical line, which is evidently positive or —* according 
as BC is < or > AC. 


It appears, therefore, in both cases, that addition is an operation ie : 


performed at B, the second extremity of AB. 

On Subtraction. 

From considerations such. as the foregoing, we shall readily con- 
vince ourselves, that to subtract one positive line from another, or 


to perform the operation a — (+ 6), can only be accomplished by 


drawing a line AC = 6 from A, the origin of the minuend, and — 


- that the equivalent difference AB — AC or a-—6 is the distance 
CB measured from the final point in the subtrahend to that of the 
minuend. 

This difference will obviously be negative wien AC is > AB, 
being always expressed by CB. | | 

_ If AC be drawn from A negatively, then CB is” 


and is , properly represented by a line equal to the sum of the two 
parts a and 84, and positive in its description. | 

It may be remarked, that this is an independent proof of the 
rule for the concurrence of signs, and one whose place is supplied 
by no other. For that rule, as proved in Arithmetic, is confined 
to the concurrence of a double operation, and in Algebra it is an 


only one for the case of the concurrence of the signs of operation 
with those of affection. 

With respect to the expressions a+(—6b) anda — (+5), itis 
evident that we obtain for them two geometrical lines, specifically 
different, but always identical in magnitude and direction; and as 


arbitrary assumption; but the proof we have here obtained is the © 


: magnitude and direction are the only elements which are denoted © 


by our algebraical system of representation, these two lines cannot 
be algebraically distinguished. 


The rules here laid down for Addition and Subtraction, verify | 


the algebraical results no less beautifully when the quantities in- 
volved are apne symbols, and the direction of the lines rect- 
angular. 

If BC | to AB be drawn upwards or dow nwards through B, 
_we shall have AC as before 
=a+(+bv—1); 
and if AC be similarly drawn through A we shall have 

CB=a—(+bV—1). 

So also, if AB be drawn in a direction inclined to the axis of a, 
and equal to a + 56“ —1, and if we draw through B, (and 


é 
| 
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therefore add) a line 
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then AC, the sum of these two, will fests bis seen to represent 


and if AB and AC represent 


atbV—1 and c+ 


respectively, we get 


CB=a—e + (b—d)¥—1. 


VIL—NOTES ON THE UNDULATORY THEORY OF 
LIGHT.—No. I]. 


THE theorems of Fresnel: respecting the motion of light in crys- 
tallized media, were demonstrated by him by means of two surfaces, 
the surface of elasticity and the ellipsoid of construction, which is 
related to the former by certain reciprocal relations. He was thus 
not only furnished with simple demonstrations of many important 
theorems, but was enabled to make use of the known properties of 
these surfaces for the purpose of suggesting corresponding sl 
ties of the rays and waves of light. 

In the present state of the theory, it may be useful to show how 
these results may be obtained with as much or more simplicity. 
from direct analysis ; and without discarding the consideration of 
these surfaces, to use them rather to give clearness and connection 


‘to our results than as methods of demonstration. 


Since the publication of Fresnel’s Memoir in the 7th volume of 
the Memoirs of the Institute, the chiet addition which has been made 
to his results is the beautiful discovery of conical refraction by Sir 
W. Hamilton: this will be shown to flow naturally from the equa- 


_ tions we get in investigating the equation to the wave surface. 


The method of proving the existence of the circle of plane con- 
tact on the wave surface is due to Mr. Greatheed: the expressions 
for /, m and m in terms of v,, v, are given by Professor Sylvester 
in the Philosophical Magazine. 

As in the preceding part, we shall take as coordinate axes the 
three axes of elasticity, and we shall assume that the forces excited 
by displacements equal to unity in those directions are a’, 6’, c® 
respectively. We shall call a line or direction, the cosines of the 
inclinations of which to the axes are /, m, n, the fine or direction 


m,n]. 
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ae displacement, is also perpendicular to the front of the wave, it may, — | 


other, 


CS h), or 


(1) Let then lx+-my +nz=0 be the equation to the plane wave, 
and let (a, , y) be the direction of vibration ; a displacement equal 


to unity will give rise to forces a*a, b*/, c*y parallel to the axes. 


Resolve these forces along the displacement and perpendicular 


to it, and let the latter direction be determined by the cosines 


a, y, so that (a, (a’, y') being perpendicular to each 


a’2 4. 2 —]1. .. (2). 

To determine a’, we arising from 
this, that the resultant, the displacement, and the line (a’, (, y') 
being all in the same plane, must all be perpendicular to some Tine 


fa + + hy =0 
fal + 9B + hy =O (A). 
fara + + hey = 0 


Eliminating h by cross multiplication, 


If (a’, be perpendicular to plane, 
and making these substitutions, the equations (1), (2), (3) will de- — 
termine a, y- 
The part of. the resultant resolved along the displacement i is 
ara? 4 4 


and, as the other part, which is perpendicular to the direction of 


as we have shown before, be neglected, we have a wave proceeding 


with a velocity 


= Varad + 4 ey? (4). 
ry If we multiply together the equations (1) and (3), put under 
the form 
| la + m3 = =— ny | 
+ (c? — a?) = — (a? B2) (B), 
| 


we the equation 


a —c 
a quadratic in =, the last term of which, ——_, is negative, and 


therefore the roots are real; the same is the case in the equation we 


should find for 4: and thus the above equations give two real di- 


rections (a, 


4 .* 
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Let Bos Yor be the roots 


62 — ¢ 
[3,35 
b2 — a? — 6?’ 
hence aa, + = O 


and the two directions are perpendicular to each other. 


_ Whatever then be the nature of the vibration of the particles in 
the plane wave, we may resolve it in these two directions, and we 
shall have two waves, whose vibrations are perpendicular to each 
other, and proceed with different velocities v, and v,, determined 
from equation (4): that is to say, any plane disturbance will pro-— 
duce two plane waves polarized in planes at right angles to “_ 
other, and of different refrangibilities. 


(3) v, v, are the roots of a quadratic equation, which may 
easily be found from the equations (A) by substituting /, m, m for 
a’, 3, y- The method is the same as that in page 7 of our * last } 
7 Number, and the result is the same, namely, 

[2 | 


at 


(4) Having determined the velocities of propagation of a plane 
disturbance, the form of the wave surface, or that surface into 
which a disturbance at any point will diverge, may be found. 


If we suppose the disturbance to be caused by an infinite number 
of plane waves passing through the origin, each too faint to cause 
any impression on the senses of itself ; then each of these waves will 
be propagated with the velocities we have just determined, and the 
form of the sensible disturbance at any instant will be the locus of 
their ultimate intersections, or the surface touched by all the ages 


‘Waves. 


At the end of a unit of time the equation to one of these planes 
1s 
| lx + my + = 
with the relations | 


The process of elimination will be found in the Cone 
Transactions, Vol. vi. Part 1. 


(5) The following three equations are found after eliminating 
the differentials, 


4 
C= 
3... 
‘ 
* 
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| 
mv 
(b? — r*) >. --(D), 
nv 


in which r= = + 2’, 


and the resulting equation to the wave surface 1 is 


If X, Y, Z be the angles which a radius vector makes with the 
coordinate axes, the polar equation is | 
a* cos? X , 6? cos? Y , _ 


6. We shall now distuss some of the equations we have ob- 


tained. 


_ The two equations (B) will in general determine two of the 


quantities 4% ”, and the directions of the vibration lines will be 


definite ; but if one of the quantities J, m, n, as for instance m, is 


zero, the two equations will become identical, provided 
| n2 (a? 6) (2 (6? — = — (10), 
and thus the direction. of the lines of vibration become indetermi- | 


nate, if 
b? — c? 


These values of J and n are real, provided a, b, ¢ are in the order . 
of magnitude, but not otherwise; and thus we should not get a 
similar property by making /=0 or n= 0. 


Putting a? — + — v* and — v* — 
and c* — v® in equation (7), we get 
— 4 — v?) x 
[n? (a? — — /? — c?) + m? Sa? — — — 
—- m? (a? — 6?) — c?) = 0 ()), 


(6? — c*) for a? v’? 


which on our supposition reduces to 


(0? v2)? — — 

Thus the two values of v? become equal to each other and to 6”, 

There is no separation of the plane wave; and therefore for the 
two directions determined by equations (E), and for these only, a 


single plane wave may be propagated. The normal to this plane is 
then a direction or line of single normal velocity. 


— 

& 

¢ 

— 

y 
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7. If sages (9) be put under the form 
cos? X cos’ Y cos* 
it becomes exactly similar in form to equation (7), and we at once 
infer that the two values of r? become equal to each other and to” 


provided 
cos ¥ = 0, Z. = 


and the radius vector, which in general cuts the ae in two 
points on the same side of the origin, will in this direction cut it 


only in one. This then is a direction or line of single ray 
velocity. 


cos X = 


8. These two directions, to which we have given Sir William — 
-Hamilton’s names, the lines of single normal velocity, and the lines 
of single ray velocity, are nearly coincident in most crystals. They 
were not distinguished by observation before Fresnel’s theory, but, 
from some remarkable optical. properties which they possess, they 
received the common name of optic axes. 

When a, 6 and ¢ are different, the. equations (E) and (F) give 
two lines inclined at equal angles to the plane of xy, and crystals 
of this kind were called biaxal ; if b= a or =e, the two directions 
coincide, and such crystals were called uniaxal. Fresnel has given 
the name of optic axes sometimes to one set and sometimes to the 
other, and considerable confusion has arisen from the ambiguity. 
We shall use the terms wave axis and ray axis; names which pre- 


vent ambiguity, and which appear ary to explain their own 
meaning. 


9. By comparing equations (7) and (E) with equations (9) and 

' (F), we see that, supposing the direction of a ray to be the same as 

that of a wave normal, the equations for determining the wave 

velocity v and the directions of the wave axes, are changed into 
those for determining the ray velocity 7 and the directions of the _ 


ray axes, by substituting a for v and for a, 5, ce, their reciprocals. 


And hence, that every relation between the wave velocity and the | 
directions of the wave normal and wave axes, will give a corres- 
ponding relation for the ray. 


If we have 
» (1, m, n, v, a, b, ¢, direction of wave axis) = 0, 


4 
& 
M 
- 
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we derive at once the _— relation 


] 
X, cos cos Z, direction of ray axis) 


10. Equation (11) gives us the inclination of the wave normal 
to the axes in terms of the two wave velocities V1» Up} for as the 
last term is the product of the two roots of the equation } 


— v,*) — v >) =—m (a? — b?) (b? — c?), 
= 


(a?—b?) (a2 — =// ce?) ..(H1). 


_ From these on, we may find 1, ’, the angles between the 
wave normal and the wave axes, in terms of the velocities v,, v,, 


Substituting the above values of / and n, this becomes 
squaring, putting a? — e* — (v,* — c*) for a? — v,2, and making a 
similar substitution for v,? — c?, the equation reduces to: 


(a? — c?)2 cos? (a*— 
(a? — c*) sine = V (a? — 0,2) (v,2 — e2) — — 0,2) (v2 — c?). 


The expression for sin ¢’ will only differ by the sign of one of the 


radicals; the product of the two expressions will therefore be the 


difference of the squares of the radicals, or 

(a2 — sine sin’ =. + (a? — — 

or — v2 = + (a? — co?) sine sine... .....(12). 

If «, « be the angles between the ray and the ray axis, we deduce © 
from the above at once | 
(-, 5) sin coe cee (13), 
or the difference of the reciprocals of the squares of the ray velo- 
cities is proportional to the product of the series of the inclinations of 
the ray to the ray axes. ‘This law was discovered experimentally by 
Biot. In Airy’s tract on the Undulatory Theory, No. 123, this 
property is inadvertently stated as belonging to the inclination of 
the ray to the wave azes. 


11. We shall now proceed to consider the form of the y wave 
surface given by er (8). 


| 


‘ 
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_ The equations to its intersections with the coordinate planes will 
| consist. in each case, of a circle and ellipse. In the plane of avy the - 
_ circle is contained within the ellipse; in the plane of yz the ellipse is 
contained within the circle; in the plane of xz the ellipse and circle 


~ eut each other, — given by the equation 


~ 1) = 0. 


The.point of intersection of the two curves is the extremity of the 
ray axis. The point where a common tangent to the two curves: 
touches the circle is the extremity of the wave axis. | | 


Fig. (2)i is a section of the surface by the plane of xz, D is the 
extremity of the ray axis, and F, the point where the common tan- 
gent EF touches the circle, the extremity of the wave axis. 


12. The equations (D) determine the point at which a tangent 
plane of given position touches the wave surface; or if the point of 
contact be given, they determine the position of the tangent plane: 
that is, if 4, m,n” and v be given, they determine 2, y, Z; con- 
versely, if x, y, z be given, the equations @) with 


la + my =», 
will determine J, m, n, v. 

If for any values of J, m, n, v, one of these equations is ‘tatishied 
independently of the values of 2, y, and z, these equations will de- 
termine, not a point, but a curve of contact with the tangent plane. 
_ Also, if for any values of 2, y, z, one of the equations is satisfied 
_ independently of /, m, n, v, we shall get only two relations between 
1, m, n, and there will be an infinite number of tangent planes at 
that point, giving a cone of contact. 


If m= = bY, the second is satisfied inde- 
pendently of 2, y, z: this is the case we have before considered of 
the plane of single normal velocity. 


The equations (D) are the equations of three spheres; hence, any 
‘two of them determine a circle: and thus it appears, that that plane - 
wave which is propagated without separation, or along the wave 
axis, touches the wave surface in a circle, the — and magni- 
- tude of which we proceed to find. 7 


We have, as before, 


substituting these values, the two equations (D) become 

* = — 6) (a? — &) 

| 

‘= ae ay 


t 
- 


84 Notes on the Undulatory Theory of Ligh —No. I. 
or putting a! +y + 2 for 7’, 


ga 4 — 62) (a? — c?) 0, 


If. we eliminate 2, We get e the projection on the plane of ay 


the equation to an ellipse, of which the axis parallel to y is the 
radius of the circle of contact, : 


This — also gives one coordinate of the centre of the circle, 


2_ PR +02 


a2 — & 


Similarly we should find 


_ And thus the position and magnitude of the circle of contact is 
completely determined. 
From the double sign + which these wives involve, we see 
that there will be four such circles on the wave surface, one at each 
end of the two wave axes. 


the second equation is satisfied indepen- 
dently of J, m, nm: this is the case at the extremity of the ray 
axis, where 


Substituting these values, the first of equations (D) becomes 


2 2 
a ee Cc a — UV 


similarly, v? — c? = Via — 
and eliminating v, we get the equation | | 
12 (62 — c2) 4 m? (a? — c?) 4 n? (a® — 6?) 


| —c gg? 4 
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which, with the equation /? + m?-+ n? = 1, gives the relations 

between 1, m,n. | | 
At this point there are an infinite number of tangent planes, 

whose ultimate intersections give a tangent cone. The equation 


found above, where a, y, z are put for J, m, n, is the equation to — 
_ the surface described by a line always normal to these planes. 


15. To find the equation to the tangent cone, call the above 
equation wu = 0, lx + my + nz =0 is a plane through the 
origin parallel to a tangent plane. _ 


We have, therefore, 
+- ydm + zdn = 0, 
ddl + mdm + ndn = 0, 


du du 


therefore + + pe = 0, 


| dm + vy = 0; 
du 


Multiplying by 4 m, n, and adding, X=0. We have then to 
eliminate p, 4,m,n between | 


u =O, + pr = 0, py = 0, + pz = 0: 
this is easily done, and the resulting equation | 
_(@— &)y 92: 


LZ 


is the equation to the tangent cone, the origin being at its vertex. 


This is an elliptic cone; the axes of any section made by a plane 
perpendicular to its axis, may be found by transforming the co- 
ordinates so that the term involving xz may vanish. 


The angle of the cone in the plane of wz may be found from 
either of the equations (17) or (18). If we make m = 0 in (17), 


l 
and solve for —, we find See 


n 


7 
‘ 
| 
4 
| 
> 
4 
~ 
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and the tangent of the angle sought 


_ difference of roots 
~ | + product of roots 


ac 
V (a — (0? — c*) 


At the extremity of the ray axis there is a cusp, or rather a 
nodal point, somewhat resembling the vertex of a very obtuse 
_ double cone. Perhaps a clearer conception of this surface may be 
formed by supposing an ellipse, as in fig. (3), to revolve round a 
diameter AB, which is not one of the axes. At A and B there will 
be two nodal points, each having an infinite number of tangent 
planes forming two cones of contact; and the two tangent planes 
FF’, GG’ will touch the surface in two circles. At the points D — 
and E there is a cusp or node of a different kind: here there will 
be only two different tangent planes, and the tangent planes HH’, © 
KK’ will touch the surface in two points. The form of the surface 
~ at A and B will represent the form of the node of the wave surface 
very nearly, but not exactly, unless the ellipse be supposed to 
change its form while revoiving. . 


From the existence of the four tangent cones at the extremities 
of the ray axes, and the four circles of plane contact corresponding 
to the wave of single normal velocity, Sir William Hamilton draw 
the following inferences, which have been fully confirmed by the 
experiments of Professor Lloyd. 


1. That if a ray be refracted into thie onal so as to proceed 
along the ray axis, then, as the position of the tangent plane at the 
extremity of the axis is indeterminate, it will at emergence be re- 
fracted into an infinite number forming a cone, and giving external © 

conical refraction. 


2. That if a single ray be incident externally at an angle corres- 
ponding to the wave axis, it will be refracted internally into an 
infinite number forming a cone, and giving internal conical refrac-— 

tion. If the faces of the crystal be parallel, the rays will emerge 
parallel to the incident ray, oeming & an external cylinder of rays. 


A. S. 


> 
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Tue name of Transversals was given by Carnot to lines considered 
_in their relations of mutual intersection. Many of their properties 
are very curious, and form interesting problems in Analytical 
Geometry, though it was not in this way that Carnot considered 
them. His method was, to proceed step by step from the more 
simple properties to the more complicated; but it seems better to 
consider each independently. | 


1. The following problem, under a lightly different form, was © 
given in one of the Problem Papers for 1836. 


If two lines AB, CD intersect in O so that AB is bisected, and 
if the lines AC, BD meet, when produced, in  E, ‘and AD, BC in 
F, then the line EF is parallel to AB. 


Take O as the origin, and AB, CD as the axes of « and. ly. 
Let OB =a, OA =-a, OD=8, OC = 


The equations to BD and AC are | 


At their intersection the equations may he combined in any man- 
ner; therefore, subtracting them, | 
y (5 ~ =) (1). 
The equations to BC and AD are 


At their intersection, subtracting them, 


- which is the same as the equation (1), and therefore is the equa- | 
tion to the line passing through the points of intersection, that is, 
to EF; and from its form it is evident that it is parallel to the axis 


of x, that is, to AB. 


. If three lines be jun from the angles of a triangle A, B, C, 
(fig 5), through one point O, and meeting the sides of the triangle 
in P, Q, R, the sides of the triangle P, Q, R will, when produced, 
meet those of A, B, C in three points, which are in the same 
straight line. 
Take L, the point of intersection of AB and PQ, as the origin, 
and these lines as the axes of x and y. 


| 
* 
- 
‘ 
. 
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Let LB = LR=a,, LA=a,, LP = LQ = 6,. 
The equations to BC and QR are 


2 2 


At their intersection we have, by subtraction, 


— the equations to AC and PR are ee 


l, 


and at their intersection we m by subtraction, 


+9 
a, 


‘Again, combining by au oe equations ¢ to BC and AC, we 
get at their point of eal 


+2) +9 = 2 


the to AP and 
| a, a by | 
and at their intersection we have, by addition, 


which is the same as (3), and therefore is the equation to COR, 
which passes through the two points of intersection. 


If in this equation we make y = 0, we have «= 24,0; , but 


in this case x = LR = Ay; therefore | 
a 
therefore equations (1) ne (2) coincide, and represent the line 
passing through M and N, and as from the form of the equation it 
evidently. passes through the origin, the three points L, M, N are in 


one straight line. 


In somewhat a similar manner it might be shown, that if three 
circles be touched, two and two, by pairs of tangents, the points of 
intersection of these tangents are in one straight line. But we shall 
pass on to another problem. 


3. If from any point A (fig. 6) in a ve of indefinite length 
equal distances AQ,, Q,Q,, &c. be measured off, and other 
equal distances AP,, P P,, &c. in the opposite direction; and 
if Q,. Q,, Q,, &e. be joined by straight lines with a point O,, 


+ 
| 
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and Pj, P,, P,, &c. with a point O,,. O, and O, being situate in 

a line parallel to the given line; then the corresponding points ‘of 

intersection of these lines will lie in straight lines forming two bun- | 

dles of lines converging to two points O’, O”; situate in O,O,. | 
Take A as the origin, AO, as the axis of z, AQ as the axis 

of | 

‘Let AO, = m, 00, =¢, AQ=a, AP = 5, 


- Let one of the fines O, P,, meet the axis of z ina point a,; then 
by similar triangles it is easy to show that 


A,  m 


The equation to 0,Q, is 


Combining them by subtraction, after multiplying by 7 and s, 


mb 


_ which gives a relation between the coordinates of the point of in- 
tersection of the lines. We shall arrive at the same result for all 
pairs of lines for which s—r is the same; consequently equation (3) 
will represent a line passing through the point of intersection of 
each pair. NOW; let x =m, then 

| ahich being independent of r and s, will be the same for all the 
lines represented by (3) for which s—~, is different, and there- 
fore they all converge to one point, which it is easy to see lies 
between O, and O,, as the value of y is less than that of ec. If we 
make = 0 in (3), we get | 


and as s—~7 increases by unity, the distance between the points 
where the different lines diverging from O' cut the axis of y is half 
the harmonic mean between a and 6. | 


Tf, instead of subtracting, we add equations (1) and (2),9 we get 


(7 +s)— + +S 


a 


which, as before, may be shown to be the equation to a line passing 
through the intersections of all lines for which r+4s is the same. 


m Ta 
and the equation to O.P,is | 
ON 
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And if r+s vary, we find, as before, that the different lines 
converge to a point 0’, such that 
= — 


b—a 

If b> a, O" is situate in O,O, produced. If a= b, O,O" 
becomes infinite, and the lines are then parallel to AQ. If b <a, 


_ the point O" is above O;. The distance between the points in 


which the lines diverging O” cut the axis of y is 


b 
It is easy also to show, that all the lines measuring from the points 
of convergence are divided harmonically : the equation to any 
line, as 02Q;; is | 


CXL 


| m * msb- sa 
| 
when y 1+ <) 1; 
similar! 


and so on for the others; and as these are in arithmetic progression, 
O,s, O,7, &c. are in harmonic progression. Instead of taking O,A 
as axis, we might have taken. any other line, as O,Q,, and the 
result would be the same. And in like manner we might proceed, 
taking O, as the origin, as also O', O", and, as the cases are similar, 
we should get the same result. Hence. the property holds for all 
the lines in the figure. For other problems in Transversals the 
reader is referred to Carnot’s Memoir, or to a small work by 
Brianchon on the same subject, where he will find some practical 
‘applications of the theory. | 


F. 


1X._ON A NEW FORM OF EQUILIBRIUM OF 
A REVOLVING FLUID. 


In a letter addressed to the French Institute in 1834, M. Jacobi 
announced the curious discovery, that the form of equilibrium of a 
fluid mass, consisting of particles attracting each other according to 
the law of the inverse square of the distance, and revolving about 
an axis, may be an ellipsoid with three different axes, the shortest 
being the axis of revolution. A demonstration of this theorem is 


4 
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given by M. Liouville in the Journal de U Ecole Pain, 
vol. xiv.; but as this work is not easy to be procured, we think a 
notice of this remarkable discovery will be acceptable to many of 
our readers. The following demonstration differs considerably from 
that of Liouville. 


Let the fluid mass revolve about the axis of z with an angular 
velocity w, and let X, Y, Z be the resolved attractions on a + particle 
at (a, y, 2). 

The equation to the free esitline of the fluidis 

(X — w’x) dx + (Y — w*y) dy + Zdz = 0. 
___ As we know that in the case of an ellipsoid X is of the form Az, — 
Y = By, Z = Cz, where A, B, C are independent of the coordi- 
nates of the point attracted, we see that this is in the form of the 
differential equation to an ellipsoid whose axes are a, 6, ¢, or 


d zdz | 


provided 


It remains to find whether these relations will give real values for 
the axes and angular velocity when a and 6 are different. 


_3M 


Now A = —. 
a2 + (0? — a”) 3a? 4 (c? — a*) 


Transform this definite integral into one whose limits are O and 4. 


by putting s = = , and we get 
du 
(a? 


calling the quantity under the radical D, 


similarly B + ud’ C= ra 
Eliminating w* and “ we get the relation 


a*b? (B — A) — — 6°) =0,. 


» 
a? 
| 
* 
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If abe different from 3, the relation between the axes must 
satisfy the equation 


0 : 
Supposing a and 6 given, this may be considered as an equation 


to determine ¢, and as it is negative when e=0, and positive when 


c= = there will be at least. one real root os Cc. 


As the above integral is zero, and 7, a is always power’ and 


] u 


ot op 


is positive if « be sina enough; when wz is small, this factor must 
be negative, or | 
a? 
therefore c? is less than a® or 6°. 
To determine the angular velocity, we have 


(a* — Aad — Bo 


he / udu 
= 5M 
0 


if a be different from 6, 


= Merge 


a, positive w is a possible quantity, and there 
being at least one possible value of c, the problem is possible. 


H. M. 


X.—RADIUS OF ABSOLUTE CURVATURE. 


Tue following is a simple method of finding the expression for the 
radius of absolute curvature, taking the arc as the independent 
variable. Let PQ, QR (Fig. 4) be two consecutive elements of 
the curve, and let us suppose them to be equal; this is the same 
thing as making ds constant, or s the independent variable. If we 
complete the parallelogram PQRS, the centre of curvature will lie in 
QS produced, which bisects the angle PQR, and if O be the centre 
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of curvature it will be the centre of a circle passing through P, Q,k | 
so that putting QO =p and QS=), © | | 

= ds? = dp. 


Let a, 3, f be the angles QO makes with the axes, then projecting 
QS on the axis of 2, we have ~— | 


a QS cos a = RT cos a = M'N’. 
Now M'N' = NN’ — NM'= MN— NM’ = 
me therefore \ cos a = dx; 


whence cos * 


Then from the relation cos? a + cos? 3 + cos? y= 1, we find 


V de” d# 


which is the well known expression. 


From these values of cos a, cos B, cos y, Cauchy’s demonstration 
of Meusnier’s Theorem follows very readily. | 

Let «=O be the equation to a surface. Differentiating twice 
with regard to 2, y, Z, we 


| 
2 
0 d?y Pz + &e. 


The terms we have left out depend Bs on the values of dx, es dz, — 

_ and therefore only on the direction in which these variations are taken, 
and will be constant if we suppose them taken along a section of 
the surface made by a plane whose intersection with the tangent 
plane is constant. If p be the radius of curvature of the section, 

Bee. a, (3, y, its s inclinations to the axes, we have 

2 

—.cosa, &e. 

| 

_ Substituting these values 


cos B du t 
ns Te (cos « + dy + cosy 7) + const 


| dx 
the cosines of the angles which a normal to the surface makes with 
du du du 
the axes are proportional to - rip ae Hence, calling 6 the angle 


= const. 


| co 
between the normal and the radius of curvature, 


When 6 = 0, p attains its maximum value = R. 


‘ 
+h 
4 
* 
: 
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o = KR cos 6, 


or the radius of curvature of an oblique section is the projection of 
the radius of curvature of the normal section which has the same | 
intersection with the tangent plane. 


H. M. 


MATHEMATICAL NOTES,. 


l. Sinci the paldiention of our first Number, we have met with 
the following proof, by Professor Jacobi of Koenigsberg, of the 
equation (K) page 16, which i is better than the one there given. 

As there, we have 


= 0 So yn) dx dy. 
Assume 


then, while wand y vary from 0 to », 7 will vary from 0 to x, 


and since z = - » 2 will vary from 0 to 1. Suppose that x and | 


yand 7, vary together, 
| dx —rdz, dy = dr. 


Substituting these values in the double integral, it becomes 


=T(m +n). (2) 


2. Diagonal of a parallelopiped. The following is an easy 
method of obtaining the expression for the diagonal of a parallelo- 
piped in terms of the edges, and the — they make with each 
other. 

Let a, b, c be the edges, r the diagonal. 

a, 2, y the angles the edges make with each other. 


A, p, v the angles the diagonal makes with the edges. 


_ Then projecting the three edges on the diagonal, we have 
r= acosd + bcosp 4 € cos». 
Again, if we project r on a, and compare it with the projection of 


the broken line a, 6, c, which joins its extreme points, on the same 
line we have the relation 


=a-+bcosy +c 


| 
+ 
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Similarly, | 
rcos p=acosy + c cosa, 

and rcos vy + bcosate. 

— by a, b, c, and add. -Then, by the first condition, 
= a? 4 b? 4 c? 4 Qab cos y 4 Zac cos 4 cosa. 


In the Ellipse a? +. 62 = b. This well-known pro- 
‘position may be very conveniently proved in the following manner, — 
The equation to a diameter conjugate to that passing through 


is 


Transposing and squaring, 
yy? 


Also, 9? = = “2 — and at the pom where the diameter 


42 
meets the curve = =, (@ Substituting these values, the 


equation becomes 


ary? (x? , 
a ae a4 
or at = 2? 4 
Similarly, P=y+y? 


therefore a? + 4 = @? re 


de 


4. pane in the papers re foe 1838. <A “iii of an ellipse is 
inscribed in a right-angled triangle having its axes coincident with 
the sides containing the right angle. Find the curve traced out by 
its centre of curvature at the point of contact with the hypothenuse, 


while the point of contact moves from one acute angle to the other. _ 


Let the origin be at the right angle, and let the sides be m, n, 


and the axes a, 6; «xy the coordinates of the point of contact, a, B 


of the centre of curvature. ti — to the ellipse is 


Since the hypothenuse i is a tangent 


yy’ 
+ = 0. 
a 
xr 
az. 
m 
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Also, we have | : 
therefore —, = 
b2 
whence — (a — ma, — = — np. 


‘Adding these, 
(= + (a? — = np. 
Again, multiplying the first by m2 and the second by n*, and 


| whence m°a + n°8 = (ma — 
the required locus, which is a parabola. 


b The area of the parallelogram formed by tangents applied — 
at the extremities of any two conjugate diameters of an ellipse 
is ‘constant. | 


Let x, y be the coordinates of the extremity of one iii: and 
a’, y' those of the other; 0, 6’ the angles which they make with a. 
Then x =a’ cos 0, y= sin 8, 2’ cos = sin @. 
The tangent whose inclination to the axis of 2 is 0, passes 
through the point a‘y': therefore, by the a in page 9, sub- 
stituting tan 0 for a, and multiplying by cos 0 


cos a sin == V(b cos 0)? (a sin 


Multiply by a’, and substitute for 2’ and y’ their values in A and 6, 
therefore | 


sin (6 — 8) = (a’ cos 4 (a@ sin 6)’ 
—= 4 ay? — ab. 
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